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ABSTRACT 

We describe a new, direct method for determining the opacity of foreground galaxies 
which does not require any a priori assumptions about the spatial distribution or the 
reddening law of the obscuring material. The method is to measure the colors and 
counts of background galaxies which can be identified through the foreground system. 
The method is calibrated, and the effects of confusion and obscuration are decoupled by 
adding various versions of a suitable deep reference frame containing only field galaxies 
with known properties into the image of the foreground galaxy, and analyzing these 
"synthetic field" images in the same way as the real images. 

We test the method on HST WFPC2 archived images of two galaxies which are quite 
different: NGC 4536 is a large Sc spiral, and NGC 3664 is a small Magellanic irregular. 
The reference frames are taken from the Hubble Deep Field. 

From the background galaxy counts, NGC 4536 shows an extinction ^4/ ~ 1 mag 
in the northwestern arm region, and lower than 0.5 mag in the corresponding interarm 
region (no correction for inclination has been attempted). However, from the galaxy 
colors, the same reddening of E{V — J) ~ 0.2 is observed in both the arm and the 
interarm regions. In the interarm region, the combination of extinction and reddening 
can be explained by a diffuse component with a Galactic reddening law {Ry ~ 3). In 
the spiral arm, however, the same diffuse, low opacity component seems to coexist with 
regions of much higher opacity. Since the exposures are shorter the results for NGC 3664 
are less clear, but also appear to be consistent with a two component distribution. 



^ Based on observations with the NASA/ESA Hubble Space Telescope obtained at the Space Telescope Science Institute, which is 

operated by the Association of Universities for Research in Astronomy, Incorporated, under NASA contract NAS5-26555. 
■^Sternwarte der Universitat Bonn, D-53121 Bonn, FRD 

''On assignment from Astrophysics Division, Space Science Department of ESA 



1 



Subject headings: galaxiesiindividual (NGC 4536) - galaxies:individual (NGC 3664) 
- galaxies:irregular - galaxiesdSM - galaxies: photometry - galaxiesrspiral - galax- 
ies:statistics - ISM:dust, extinction - ISM:structure 



2 



1. INTRODUCTION 
1.1. Historical background 

The issue of whether galaxy disks are significantly 
opaque is a lively topic of current research and debate 
(e.g., Davies & Burstein 1995, and papers therein), 
and the conclusions of different workers often differ 
widely. Statistical evidence that most large presented 
by Valcntijn (1990) based on a diameter-selected sam- 
ple of galaxies from the ESO survey. Valentijn's result 
indicated that the bright galaxies arc substantially 
opaque right out to the 25 B mag arcscc^^ isophote. 
This conclusion has been challenged by Burstein and 
others, who initially questioned the selection effects 
inherent in the sample used by Valentijn, and have 
maintained that very little extinction is present. More 
detailed studies on individual galaxies have generally 
concluded in favor of at least some extinction (e.g., 
studies of multicolor optical and infrared images by 
Block et al. 1994, and Pclcticr et al. 1994). There 
is also some indication that even galaxy halos may 
contain dust; Zaritski (1994) compared faint galax- 
ies in fields bordering nearby spirals (NGC 2835 and 
NGC 3521) and in control fields further away on the 
sky, and found a small difference in the average col- 
ors. However, from analysis of surface photometry 
on a sample of 1700 galaxies, Giovanelli et al. (1995) 
have recently concluded that less luminous galaxies 
are completely transparent, and that even the most 
luminous spirals suffer from only small amounts of ex- 
tinction (0.1 mag) at an isophotal radius of 23.5 mag 
arcsec"^ in /. 

Tests of the opacity of spiral disks which rely on 
intrinsic properties of the galaxies (surface brightness, 
inclination angle, isophotal profiles) are strongly de- 
pendent on the necessary assumption of just how the 
dust is distributed within the galaxy (Witt, Thron- 
son, & Capuano 1992; Witt & Gordon 1996). For 
instance, layered models with a physically thin but 
optically thick dust layer behave like optically thin 
disks (Disney, Davies, & Phillips 1989). In the case 
of edge-on galaxies, surface photometry can help de- 
termine the geometrical dust distribution (Kylafis & 
Bahcall 1987), but in more general cases the results 
are inconclusive. The use of stars or star clusters em- 
bedded in the galaxy can only tell us something about 
the remaining opacity between us and the star, and is 
clearly biased to low opacities. HII regions have also 
been used; they have the same bias, but in addition 
those bright enough to be easily observed may alter 



the surrounding ISM so drastically that they are poor 
test objects (e.g., Allen et al. 1997). 

A more direct approach has been proposed by 
White & Keel (1992), and followed up by White, Keel, 
& Conselice (1996) and Berlind et al. (1997). They 
have looked at partially-overlapping galaxy pairs, the 
first two studies at optical wavelengths, the last one in 
both the optical and the near-infrared. Assuming that 
the background galaxy isophotes are smooth and/or 
symmetric, they use the unobscured part of the back- 
ground galaxy to calculate the amount of background 
light that is lost in crossing the foreground disk. All 
three studies have found higher extinction in the arms 
than in the interarm regions, and extinction curves 
flatter than Galactic, consistent with a patchy dust 
distribution in the foreground galaxy. 

The use of background illuminating objects has 
the advantage that no assumptions need to be made 
about the relative distribution of dust and light in 
the foreground galaxy. Also, the overlapping-galaxy 
method in particular provides results with good spa- 
tial resolution. However, it is necessary to assume 
that the light distribution in the background galaxy 
image is smooth and centrally symmetric. In addi- 
tion, there is an uncertain contribution of scattered 
light from the background galaxy into the line of sight. 
But perhaps the biggest disadvantage is that there 
are only about two dozen galaxy pairs which offer the 
right geometry for the method to be applied (White 
et al. 1996). 

1.2. Distant field galcLxies as probes 

A potentially more general way to establish the 
opacity of a galaxy is to look for effects on the counts, 
colors, and morphologies of distant field galaxies seen 
through the disk of the foreground galaxy. This con- 
ceptually simple approach unfortunately suffers from 
an equally simple problem; the general distribution of 
disk stars in the foreground galaxy reduces the con- 
trast on background galaxy images, and foreground 
star clusters and HII regions add strongly to the con- 
fusion. High sensitivity is required in order to reach 
a surface density of background galaxies which is 
large enough to derive a statistically significant re- 
sult. High angular resolution and color information 
are required to assist with the separation and identi- 
fication of background galaxy images from objects in 
the foreground galaxy. 

Hubble Space Telescope (HST) Wide-Field Plan- 
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etary Camera 2 (WFPC2) images offer a significant 
advantage over images from ground-based telescopes 
in tlieir greatly improved resolution. This reduces 
the perturbation of faint galaxy images by nearby 
foreground objects, and facilitates the separation of 
stellar objects from galaxies. Color information is al- 
so important; a foreground compact HII region will 
be blue compared to the image of a distant galaxy, 
whether or not there is any amount of obscuration. 
Once galaxies have been identified, their counts and 
colors can be compared with counts and colors of 
galaxies in reference fields not covered by a foreground 
galaxy. 

1.3. The data analysis problem 

The HST WFPC2 Archive contains many images 
of background galaxy fields and of fields covering 
parts of nearby galaxies. As an example of the data 
available for this study, and an illustration of the anal- 
ysis problems which must be overcome before quan- 
titative results can be obtained, we show in the left 
panel of Figure 1 one WF frame (WF3, 1^3 x 1^3 in 
size) from a typical reference field (the Hubble Deep 
Field)|^, and in the right panel a similar-sized portion 
of a deep exposure on NGC 4536. We shall describe 
both of these images in more detail in following sec- 
tions of this paper. For the moment it is sufficient 
to note that the differences in the counts and colors 
of background galaxies visible in the two panels con- 
tain the information we seek on the opacity of the 
foreground galaxy. Only a few background galaxies 
can be discerned visually with any certainty in the 
right panel. However, before we can conclude that the 
foreground galaxy is quite opaque, we obviously need 
to understand quantitatively how the data analysis 
methods we use to determine the counts and colors 
of the background galaxies are affected by structural 
features in the foreground galaxy. 

1.4. The "Synthetic Field" Method 

In this paper we describe a new method, the Syn- 
thetic Field Method^ which permits us to quantify 
and calibrate the effects of crowding and confusion 
and to determine the accuracy with which statements 

*Noise has been added in order to approximate the shorter in- 
tegration times used for the exposures of NGC 4536, as well as 
the higher sky background at the position of the spiral. The 
mosaic is displayed in the same scale as the NGC 4536 color 
image. 



can be made about the total opacity and reddening 
of an average line of sight through the foreground 
galaxy. These synthetic fields are obtained by adding 
the HDF reference field images into the foreground 
galaxy image; the combined image is then processed 
in the same way as the original foreground galaxy 
image alone. The procedure is repeated with differ- 
ent values of extinction applied to the reference field, 
and a plot drawn of the final number counts in each 
simulation versus the assumed extinction in order to 
ascertain the "best" fit to the actual counts of re- 
al background galaxies, can be estimated separately 
from the biases owing to confusion and crowding in 
the foreground galaxy field. Finally, a comparison be- 
tween the average colors of, respectively, the control 
HDF galaxies and the real background galaxies will 
allow an estimate of the degree of reddening. 

Compared to the foreground galaxy fields we an- 
alyze, the HDF frames are virtually noiseless; syn- 
thetic fields can therefore be constructed by adding 
suitably attenuated versions of them directly into the 
foreground galaxy images, without degrading the sen- 
sitivity of the data. 

Figure ^ shows 4 cumulative histograms of galaxy 
counts in Ifruw images which illustrate the steps in 
the Synthetic Field Method: Histogram a is for field 
galaxies in a version of the HDF to which noise has 
been added in order to approximate the shorter inte- 
gration time used for the exposures of NGC 4536, as 
well as the higher sky background at the position of 
the spiral; 6 is for HDF "control" galaxies in an un- 
extincted NGC 4536 simulation^; c is for HDF galax- 
ies in a simulation with ~ 0.6 mag of extinction at 
and finally d is for real background galaxies in the 
NGC 4536 field. From the figure, we would conclude 
that there is about 1 mag of extinction at / averaged 
over the WF2 and WF3 fields in this region of the 
foreground galaxyj^], and about 0.1 mag of average ex- 
tinction at / in the region of NGC 4536 covered by 
the WF4 field. 



^That is to say a WF field of NGC 4536 into which the WF 

fields from the HDF have been added. 
^The / control field has been attenuated to mimic 0.6 mag of 

extinction. 

^The results for each separate field do not differ significantly; 
the numbers of galaxies have been normalized to the area of 
one WF field. 
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1.5. Limitations of the method 

Clearly, one major limitation of our method is that 
the results are only statistical in nature, i.e., the 
method does not give an answer for the extinction and 
reddening along any specific line of sight through the 
foreground galaxy, but only a statistical estimate of 
these quantities averaged over some substantial frac- 
tion of the foreground galaxy image. The accuracy of 
the method is limited not only by the counting statis- 
tics of background galaxies, but also by the intrinsic 
inhomogeneities in their spatial distribution. These 
uncertainties must also be quantified in order to pro- 
vide a convincing result. In addition, there are several 
smaller selection effects which plague the analysis and 
which will be discussed later; for instance, the fact 
that field galaxies also show a mild color-magnitude 
dependence (Tyson 1988; Williams et al. 1996) must 
also be taken into account. 

Another limitation is that, while the numbers of 
galaxies probe the whole disk (there is useful informa- 
tion in both detecting or not detecting a background 
galaxy), the conclusions from the colors of the galax- 
ies will be biased in favor of regions of lower opacity 
(because we need to detect the galaxies in order to 
measure there colors). Therefore, it is possible to de- 
termine an extinction law by combining the number- 
count and color information only when we can safely 
assume that nowhere in the region of interest is the 
extinction so high that we cannot sample it with back- 
ground galaxies (see Appendix). 

1.6. Plan for this paper 

In this paper we present a description of the Syn- 
thetic Field Method through a detailed application 
to deep WFPC2 images of the Sc spiral NGC 4536. 
In order to examine the range of applicability of the 
method, we also analyze shallow images of the SBm 
Magellanic dwarf irregular NGC 3664 (DD095). We 
find substantial opacity in both galaxies, although the 
precision of the result is low for NGC 3664, a con- 
sequence of the shorter exposure time used for the 
observations of this galaxy. In future papers we in- 
tend to apply the method to images of other nearby 
galaxies drawn from the HST archives. 

We begin in with a discussion of the reference 
field and the processing necessary in order to render 
it useful for our present purposes. Next in ^ we de- 
scribe the data on the foregound galaxy NGC 4536 
which we have chosen to illustrate the method. Sec- 



tion y describes the methods we use to identify and 
measure the magnitudes of field galaxies. Section || 
presents the synthetic field simulations and the re- 
sults on extinction through the disk of NGC 4536. 
Section ^ presents an abbreviated discussion apply- 
ing the method to NGC 3664. In §^ we examine the 
various sources of errors and biases we have been able 
to identify. Finally, we summarize and discuss our re- 
sults on the extinction and reddening through these 
two galaxies in §||. 

2. THE REFERENCE FIELD 

The choice of the reference field is clearly of prime 
importance to this method. For simplicity, images of 
this field ought to be available with the same camera 
and range of colors as the images of the foreground 
galaxy, and ought to be entirely representative of the 
sky in areas free of large nearby foreground galaxies. 
There are by now many fields observed with the WF- 
PC2 which could satisfy these criteria, including the 
Medium Deep Survey fields (typically exposures of 
> 5000 s in the F606W and F814W passbands; Grif- 
fiths et al. 1994; Ratnatunga, Griffiths, & Ostran- 
der 1998) and the field around the weak radio-galaxy 
53W002 (24 orbits in 3 colors; e.g., Pascarelle et al. 
1996a,b). However, the Hubble Deep Field (HDF; 150 
orbits in 4 colors, Williams et al. 1996) offers the 
deepest set of images. Its properties are also present- 
ly the object of intense scrutiny by many workers, 
and detailed analyses of the counts, colors, and red- 
shifts of the galaxies visible in it are becoming avail- 
able (e.g., Cohen et al. 1996; Gwyn & Hartwick 1996; 
Villumsen, Freudling, & Da Costa 1997; CoUey et al. 
1996; Madau et al. 1996). The question as to what 
extent our chosen reference field is "representative" 
can therefore be answered with considerable certain- 
ty. For example, at TOF8i4iy,AB 26 (appropriate 
to the present work), the number-counts of the HDF 
galaxies are known to be complete, and the galaxy 
number-magnitude relation is similar to that of oth- 
er fields (Wilhams et al. 1996). The HDF therefore 
provides an excellent reference field for our purposes. 

2.1. Processing of the HDF data 

We used the Version 2 Release WFC images of 
the HDF (Williams et al. 1996) retrieved from the 
STScI/HST archive. This release of the HDF was 
made from dithered exposures using the "drizzling" 
method (Fruchter & Hook 1997). There are 3 con- 
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tiguous WFC fields in this data set (and one PC 
field which we do not use here). Each WFC field 
was recorded in the 4 photometric bands designated 
F300W, F450W, F606W, and F814W (cf. the HST 
Data Handbook) ; we have not used the images in the 
UV band F336W. 

Several processing steps were carried out on the 
HDF data set in order to render it compatible with 
the images of NGC 4536 and NGC 3664. First, the 
HDF images were all re-binned to their original pixel 
size of 0"! since this is the scale of the (un-dithered) 
exposures of NGC 4536 and NGC 3664 (and most oth- 
er exposures of nearby galaxies in the HST archives). 

Second, since there is no HDF image in the F555W 
passband used for the foreground galaxies, we creat- 
ed one by interpolating linearly between the existing 
F450W and F606W images, taking into account the 
differing responses of the three filters. This proce- 
dure will produce the right result for galaxies that 
have a linear spectrum between 4520 A and 5935 A, 
and an error in the fiux at 5400 A proportional to 
the deviation of a galaxy spectrum from linear in this 
wavelength range. For redshifts lower than z = 0.13 
(v - 40000 km s"^), the 4000 A break is blueward 
of the interpolation range, and the spectra should 
be fairly linear. However, a significant fraction of 
the HDF galaxies recovered in the simulations have 
0.13 < z < 1. Most of the recovered HDF galax- 
ies have an observed (isophotal) mpsuw brighter 
than 25 (Figure ^ . From the estimates of Villumsen 
et al. (1997) for R magnitude limits, and assuming 
(i?- J) - 0.5 (Small et al. 1995), the recovered HDF 
galaxies should be at redshifts z < 1. These galax- 
ies are also mostly brighter than m.FsssH'.intorpoiatod 
= 26 (Figure ^). The majority of HDF galaxies 
brighter than mpQQQw = 26 are bluer than 1.5 in 
(i3_F450H' — ^F606Vk) (Williams et al. 1996, their Figure 
32). These estimates allow us to model the error in 
our interpolated F555W image by using stellar pop- 
ulation synthesis models (Babul & Ferguson 1996) in 
order to produce spectra of galaxies with different star 
formation rates. From the results for the synthetic 
spectra, galaxies with z < I and (i?F450W — Vpeoew) 
< 1.5 could suffer from a systematic shift of ~ -0.03 
mag in TOF555W,interpoiated- The crror will be in the 
sense of very slightly overestimating the reddening of 
the real galaxies relative to the simulations. However, 

the dispersion in (mF555Ty,intorpolatod - WFSSSVK.real) 

is ^ 0.08 mag, of the same order as the dispersion 
of the mean (mF555W,intorpoiatod - mpsnw) color be- 



tween the 3 HDF WFC fields. We conclude that the 
error introduced into the interpolated F555W image 
is negligible. 

2.2. Making synthetic fields with the HDF 

The HDF images at the full sensitivity of the orig- 
inal observations are added directly to the images of 
the foreground galaxies in order to create the syn- 
thetic fields. The signal-to-noise ratio of the HDF is 
so high that this simple addition results in an image 
which retains the noise level of the foreground galaxy 
image. Before the addition, the HDF can be attenu- 
ated and reddened by any amount in order to create 
a set of synthetic images in which one can search for 
and measure "real" background galaxies as well as the 
HDF "control" galaxies. 

In order to provide a visual demonstration of how 
many galaxies in a "representative" field are lost to 
both crowding and extinction, we have also created 
a special version of the HDF, simply by adding an 
appropriate amount of noise to it. We refer to this 
as the "degraded" HDF. After adding the noise, its 
signal-to-noise ratio is the same as that of the images 
of the foreground galaxy, which have a shorter expo- 
sure time and a higher sky background. The WF3 
frame from the "degraded" HDF image is shown in 
Figure 1 (left panel). We produced this color image 
by combining together the yF555VK,intcrpoiatod frame as 
blue, the Ipsuw frame as red, and an average of the 
two as green. The image is displayed in a square root 
scale which favors the dynamic range of the NGC 4536 
color image; this latter image is shown in the right 
panel of Figure 1. The "degraded" HDF represents 
the number of galaxies that would be recovered in 
the total absence of any foreground galaxy, and hence 
without confusion with globular clusters, HII regions, 
and stars, but at the same position in the sky and 
with the same exposure time as the images of NGC 
4536. 

The differences in the numbers of background gal- 
axies which can be recognized between the two frames 
in Figure 1 is largely due to crowding and confu- 
sion, as will be seen in § |5.1| from the comparison 
between the "degraded" HDF and the synthetic im- 
ages without extinction]^. These data, however, also 
imply a significant amount of extinction, as we shall 
subsequently show. The power of the synthetic field 



Synthetic fields produced by adding the HDF with zero atten- 
uation to the NGC 4536 frames. 
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method resides in its ability to separate the effects 
of extinction from those of confusion and crowding; 
such a large number of galaxies is lost to crowding 
only, and it has such significant effects on the pho- 
tometry (j |7.5| ), that it would be impossible to do this 
project without accounting properly for it. 

3. DEEP IMAGES OF NGC 4536 

3.1. Observations 

The images of NGC 4536 were originally acquired 
with the Hubble Space Telescope (HST) Wide-Field 
Planetary Camera 2 (WFPC2; Holtzman et al. 1995) 
for the Distance Scale Key Project; they have already 
been pubhshed by Saha et al. (1996), who give a very 
detailed account of the observations and the pipeline 
calibration of the data. The WFPC2 field is shown in 
Figure |^, overlaid on an STScI Digitized Sky Survey 
image of the galaxy. We do not use the PC frame 
(WFl) in this work. 

NGC 4536 is a bright spiral galaxy Sbc(rs)I-II (de 
Vaucouleurs et al. 1991), with i?25 = 3'79, in the di- 
rection of Virgo, at a distance of 16.2 Mpc (Saha et 
al. 1996). This yields a spatial resolution of ~ 8 pc 
per WF pixel. The Galactic extinction towards NGC 
4536 is Ab = 0.00 (de Vaucouleurs et al. 1991). 

There are 34 individual exposures in the F555W 
passband, taken at 17 discrete epochs between 1994 
June 3 and 1994 August 9; in the F814W passband, 10 
images were taken at 5 discrete epochs, over the same 
two-month period. For each filter, two individual ex- 
posures 2000 seconds long were taken back-to-back in 
each epoch, to allow for cosmic ray (CR) removal by 
an anti-coincidence technique described in Anderson 
et al. 1995 and in Saha et al. 1996. 

3.2. Further processing 

We began working on the images as retrieved from 
the archive, i.e., bias-subtracted and flat-fielded. The 
values of "warm" pixels (identified in the calibration 
pipeline as having a moderately high dark count dur- 
ing the time of observation) were repaired by replac- 
ing the measured dark current with the interpolated 
dark value between two bracketing dark frames (Voit 
1997). We then removed the cosmic rays in each im- 
age using the "CR-split" procedure (Anderson et al. 
1995). 

Since we are not looking for time variability, we 
combined all the data in each passband into a sin- 



gle deep image. Relative shifts between the differ- 
ent epochs were found with a cross-correlation tech- 
nique. The images at all epochs were nominally tak- 
en with the same telescope roll-angle; however, giv- 
en that distinguishing between point and extended 
objects is crucial for our project, we also used this 
technique to check that in fact there were no un- 
intended relative rotations between the frames. Fi- 
nally, the images were corrected for geometric dis- 
tortion (Trauger et al. 1995), ahgned, and co-added 
using the "Drizzle" variable-pixel linear reconstruc- 
tion algorithm (Fruchter & Hook 1997). The result- 
ing images have total exposure times of, respectively, 
68 X 10^ seconds in the F555W passband, and 20 x 10^ 
seconds in the F814W passband. 

4. BACKGROUND GALAXY 
PHOTOMETRY 

4.1. Object extraction 

We used SExtractor (Bertin & Arnouts 1996) to 
perform the Ipsuw a-nd Vfsssw photometry on all 
objects, up to the 24.5 mag^^ arcsec"^ isophote in 
IpaiiW- The parameters for the search and photome- 
try were chosen to optimize the background determi- 
nation, as judged from the lack of "holes" (oversub- 
tracted regions) in the background-subtracted NGC 
4536 images. Likewise, the threshold surface bright- 
ness was selected so as not to "find" noise. We 
performed the search and fixed the threshold in the 
F814W passband not only because of the actual se- 
lection bias in favor of red objects (see below §7.1), 
but also because our photometry was limited by the 
shorter exposure time of the Ipauw frame. 

From the list of all objects found by SExtractor, we 
selected the background galaxies by visual inspection 
of the images. We used the morphology, the fuzzi- 
ness and, especially in crowded regions, the color. 
In crowded regions, we selected preferentially fuzzy 
objects with {Vf555w — Ipauw) > 0.5. Our search 
and selection criteria mean that, in principle, from 
the number-counts we will measure the extinction at 

If814W- 

Figure ^ shows the cumulative histograms of real 
galaxies found in the background of NGC 4536, in 
the F814W passband. The left panel of this figure 
displays the results for the WF2 and WF3 fields com- 
bined (the numbers are normalized to the area of one 
WF field); the right panel shows the results for the 
WF4 image. We analyzed the WF2 and WF3 images 
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of NGC 4536 together, given that they span most of 
the northwestern arm of the galaxy. In contrast, the 
WF4 field covers the northern interarm region and 
a small area that appears to be beyond the visible 
disk of the galaxy. We found a total of, respectively, 
14 and 16 real background galaxies in the WF2 and 
WF3 fields, and 35 galaxies in the WF4 field. The 
difference between the number-counts in this figure 
and the number-counts of the HDF galaxies (Figure 
^ is caused by the combined effects of systematics, 
crowding and, finally, what we are interested in: ex- 
tinction. 

5. SIMULATIONS AND RESULTS 
FOR NGC 4536 

Synthetic fields were produced with varying 
amounts of (spatially uniform) attenuation and red- 
dening (Table 2). The idea is to attenuate the HDF 
fields until one recovers the same number of simulated 
and real background galaxies. This attenuation, de- 
rived from the number-counts, should correspond to 
the average extinction through the disk of NGC 4536 
to an uncertainty given by the Poisson and clustering 
errors. In practice, what we measure is: 



Aj ^ -2.5- C-log{^). 



(1) 



A^o is the normalization, in this case the number of 
HDF galaxies recovered in the synthetic fields without 
attenuation, and C parametrizes the dependence of 
the number vs. extinction curve on object selection 
biases (§ ^.3|) . 

The reddening suffered by the background galax- 
ies is found by comparing their average color to the 
average color of the HDF galaxies recovered from 
the synthetic fields. Because galaxies have a color- 
magnitude dependence, one must be careful to per- 
form this comparison with simulations that have the 
same extinction as that measured (from the counts) 
for the real galaxies or that have the same magni- 
tude limit as the extinction-corrected limit of the real 
galaxies. The selection biases of the real and syn- 
thetic samples of galaxies, and especially the color 
selection bias (§7.1), should also be the same for the 
comparison to be meaningful. 

To produce the synthetic fields, we: (1) scaled the 
HDF WFC frames to the exposure times of the NGC 
4536 mosaics; (2) attenuated the HDF frames to mim- 
ic different amounts of attenuation and reddening; (3) 



added them in directly to the images of the spiral. To 
improve the statistics, at each attenuation/reddening 
we added the three HDF WFC frames, one at a time, 
to each one of the NGC 4536 WFC frames. What 
we call one simulation includes the results of repeat- 
ing the procedure with 3 synthetic fields, one for 
each HDF WFC frames^. As in the case of the real 
background galaxies, the photometry on the synthet- 
ic frames was done with SExtractor, and the thresh- 
old surface brightness was fixed in the / image. The 
identification of control galaxies was done, like that of 
real background galaxies, by visual inspection of each 
synthetic frame (as opposed to, for example, through 
cross-correlation with a list of known HDF galaxies). 
Systematics and confusion should thereby affect syn- 
thetic and real galaxies in the same way. 

Since we did not know a priori the extinction law 
for the disk, we made simulations with both Galactic 
reddening (Fitzpatrick 1986) and completely grey ex- 
tinction (i.e., the HDF was attenuated equally in both 
the F555W and the F814W passbands). We hypoth- 
esize that the extinction through galaxy disks will lie 
somewhere in between these two limiting cases. An- 
other reason to have the two sets of simulations is that 
they should establish a relative scale that takes into 
account all possible systematic errors, against which 
to measure the attenuation and the reddening of the 
background galaxies. We also constructed synthetic 
field simulations with a "mixed" extinction law (the 
HDF Vf555w frame was more attenuated than the 
Ipsiiw frame, but less than for Galactic reddening), 
to check that the color scale established by the two 
limiting cases was roughly linear. We covered a large 
range in attenuation, from no attenuation at all up 
to 2.6 mag in both passbands. Since NGC 4536 is in- 
clined at 65° to the line of sight (de Vaucouleurs et al. 
1991), we performed these highly-attenuated simula- 
tions anticipating the possibility of a face-on optical 
depth of TV, Johnson — 1, both with grey and with 
Galactic reddening extinction curves. Table 2 shows 
the input parameters of the simulations, as well as the 
number and average output color of the HDF galaxies 
recovered from the NGC 4536 frames; a line indicat- 
ing the actual NGC 4536 observations is included for 
comparison. 



^In all the figures, the numbers of recovered HDF galaxies are 
averages of the 3 synthetic fields. 
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5.1. Completeness 



5.2. Real vs. control galaixies 



To assess the effects of crowding alone, we com- 
pared the "degraded" HDF frames (§2.2) to the simu- 
lation without extinction, i.e., a synthetic field where 
the HDF was added to the NGC 4536 frames without 
any attenuation. Figures || and ^ compare the aver- 
age numbers, of HDF galaxies only, recovered from 
the NGC 4536 frames in the simulation, against the 
average number per mpsuw and mpsbbw bin, re- 
spectively, of galaxies in the "degraded" HDF. The 
effect of crowding on the number-counts is, indeed, 
dramatic. 

Given the isophotal Ipsuw brightness limit, while 
the HDF galaxies in the "degraded" frame are com- 
plete to at least mpsuw = 25, when added to the 
NGC 4536 they are complete only to mpsuw = 24 
in the WF4 frame and, optimistically, to mpgnw — 
23 in the WF2 and WF3 frames. This is the range 
where galaxy clustering is most noticeable (§7.4). On 
the other hand, we were able to detect galaxies up to 
an integrated (isophotal) magnitude of mi?8i4W = 26 
in all frames. The Vp^ssw histogram shows that we 
were able to detect many more blue, faint galaxies in 
the WF4 NGC 4536 field. 

The histograms of the recovered synthetic galaxies, 
however, might suffer from two competing effects: the 
confusion of faint wings with the foreground disk, and 
the contamination of the photometric aperture with 
foreground objects. As described below (|7.2|), when 
using isophotal magnitudes galaxies are measured, on 
average, 0.1 mag too bright in the F814W bandpass, 
and up to 0.6 mag in the F555W filter, depending 
both on the crowding and on the faintness of the ob- 
ject. Fortunately, our investigation does not require 
absolute or total measurements of the galaxy bright- 
nesses, as long as everything is compared in the same 
relative scale provided by the simulations. 

An absolute measurement of the total magnitude 
limit is needed only to estimate the expected cluster- 
ing error. For this purpose, we correct the mpsiiw 
completeness limit for the faulty background subtrac- 
tion, and we also assume that total magnitudes are 
0.3 brighter than isophotal. After an additional cor- 
rection for the extinction (§^.3[), we set our complete- 
ness limits at mpsuw — 23.7 in the WF4 field, and 
at mpsiiw = 21.7 in the WF2 and WF3 frames. 



Figures ^ |ri|, and |2| compare the number of re- 
al galaxies to the number of HDF "control" galaxies 
recovered in the simulations without attenuation, per 
If814W (isophotal) magnitude bin, Vf555w (isopho- 
tal) magnitude bin, and color bin, respectively. The 
effects of extinction are readily noticeable in the 
WF2 and WF3 fields, both from the reduction in the 
number-counts and from the shift to fainter magni- 
tudes of the real galaxies. For the WF4 field, any ef- 
fects are much less apparent. The reddening, though, 
is harder to discern from Figure ^2|, even for the more 
extincted WF2 and WF3 images. 

The relative photometric error for the individual 
galaxies was estimated by comparing the three mea- 
surements we had of each one of the simulated objects 
in the synthetic fields without attenuation (each HDF 
frame was added to each NGC 4536 frame). There- 
fore, it includes the effect of differential crowding, and 
amounts to 0.16 mag at Ipsuw (isophotal) and 0.26 
mag at Vf555w (isophotal). However, the error in 
(yF555W — If814:w) for each galaxy is only 0.10 mag, 
given that the color has been measured using a small, 
fixed aperture of 0"5 diameter, precisely to minimize 



the impact of crowding (see §7.2) 



5.3. The extinction measurement 

Figure |l^ shows the comparison between the num- 
ber of real galaxies (horizontal lines) and the average 
number of "control" galaxies recovered in the simula- 
tions at each given Ipsuw attenuation. The top pan- 
el displays the results for the WF2 and WF3 fields; 
the bottom panel, those for the WF4 image. The 
solid triangles are for simulations with Galactic red- 
dening; the open triangles, for synthetic fields with 
grey extinction. The different shapes of the curves 
for both types of extinction illustrate the dependence 
of the number-counts on the colors of the galaxies. 
The actual measurement of the extinction Ai was 
performed by linearly interpolating between the sim- 
ulations bracketing the real galaxies. The error bars 
include Poisson and clustering errors (§7.4). 

We find for the arm region of NGC 4536 (WF2 and 
WF3 fields) an absorption of ^/ = 0.74 ±0.49 mag in 
the case of grey extinction, and Aj = 1.07 ±0.49 mag 
in the case of a Galactic reddening law. Without the 
contribution for clustering, the purely Poisson error is 
±0.26 mag. For the interarm and outside area (WF4), 
the results are, respectively, Aj < 0.46 mag (< 0.33 
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mag without clustering error) for grey extinction, and 
Aj < 0.52 mag (< 0.39 mag without clustering error) 
for Galactic reddening. The +la error of 0.37 (0.21) 
mag is already included; the figure shows that the 
difference between the attenuation inferred from the 
two models increases with extinction. The difference 
is in the sense that we expect; it is due to our col- 
or selection bias, which favors the detection of faint 
galaxies (statistically blue) when they are reddened 
(§7.1). The difference, however, is smaller than the 
errors. 

5.4. The reddening measurement 

For the background galaxies, we measure a mean 
color of < (V — I) >ivF2&3= 1-36 ± 0.14 in the arm, 
and of < (V - I) >wf4= 1.40 ± 0.12 in the interarm 
and "outside" region. Figure |lj shows the compari- 
son of the color of the background galaxies to those of 
the simulations. We remind the reader that this is the 
appropriate comparison, rather than a confrontation 
between the two samples (i.e., WF2-hWF3 and WF4), 
given their different magnitude limits and all the sys- 
tematics that come into play in the presence of severe 
crowding. We also point out that the purpose of this 
figure is to compare the color of the background galax- 
ies to the synthetic fields closest in extinction and 
object selection biases; the locus of the real galaxies 
in the plot does not necessarily imply an extinction 
law for the obscuring material in the foreground disk, 
unless we can safely assume that the whole region 
is sampled by galaxies whose colors we can measure 
(§ ^.3| and Appendix) . The left panel shows the results 
for the WF2 and WF3 fields; the right panel, those for 
WF4. The simulations with Galactic reddening (open 
triangles) do not follow exactly the theoretical extinc- 
tion law (solid line) for large absorptions, neither do 
the simulations with grey extinction (filled triangles) 
follow a constant line in color. The departure of the 
grey simulations from the zero-r edd ening line is due to 
our selection for red galaxies 



7.1 



and to the color- 
magnitude relation of galaxies. As we go to larger ex- 
tinctions, given our detection magnitude limit we are 
left with statistically brighter galaxies, which are also 
intrinsically redder. The deviation of the reddened 
simulations from the theoretical Galactic reddening 
line is due to the systematic error in the bac kgro und 
subtraction in crowded areas described in §7.2. In 



problem. This is also the cause of the relative shift 
between the scales of the two panels. The WF2 and 
WF3 images of NGC 4536 are more crowded than 
the WF4 field; statistically, all galaxies in the former 
fields have been measured relatively bluer, producing 
the downward shift of the whole color vs. extinction 
plot. All these systematic effects will affect the real 
galaxies in the same way and should not bias the esti- 
mate of the reddening, as long as the scale defined by 
the simulations with the two limiting extinction laws 
is linear. To check that this is the case we have per- 
formed simulations with a "mixed" extinction law|^ 
The output of the "mixed" simulations is shown in 
Figure |l^. When compared with their input line, the 
simulations with a "mixed" reddening law seem more 
affected by the selection bias that pushes the grey sim- 
ulations redward than by the background subtraction 
that pushes the Galactic reddening simulations blue- 
ward. The simulations, however, follow the input line 
within the errors, which shows that the color scale is, 
indeed, roughly linear. 

The background galaxies behind the northwestern 
arm of NGC 4536 (WF2 and WF3 fields, Figure |l|, 
left) fall in between the lines for the completely grey 
extinction and the Galactic reddening law in the the 
color vs. extinction plot. Since, unfortunately, the 
color line of the grey simulations shows a kink pre- 
cisely at the location of the real background galaxies, 
to find the reddening we will rather use the color ob- 
tained by interpolating between the grey simulations 
with Aj = 0.58 mag and Aj — 1.51 mag, at the ex- 
tinctions (inferred from the grey and Galactic redden- 
ing models, respectively) of the background galaxies. 
The interpolation yields <{V - /)>= 1.22 ± 0.11 at 
Ai = 0.74 mag, and <{V - /)>= 1.25 ± 0.11 at A/ = 
1.07 mag. Therefore, the background galaxies suffer 
a reddening of ElVpssbW - Ifsuw) < 0.32 mag (grey 
extinction model) or E{Vf555w — Ifsuw) < 0.29 mag 
(Galactic reddening law). The error in the reddening, 
already included, is 0.18 mag. 

The background galaxies in the interarm region 
(WF4, Figure |lj, right), however, seem to follow the 
Galactic reddening line. Comparing their color to 
the mean color of the simulation without extinction, 
they exhibit a reddening of E{Vf555w — Ipauw) = 
0.20 ± 0.13 mag, regardless of the attenuation model 



the case of WF4, even if it is less crowded, we still see 
the effect because our magnitude limit is fainter, and 
faint galaxies are proportionally more affected by the 



^'^That is, the V light of each individual galaxy is more attenuated 
than the / light, but less than in the Galactic reddening case; 
this type of extinction can be caused by an unresolved clumpy 
dust distribution (Witt & Gordon 1996) 
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used. 



magnitudes|^. 



6. NGC 3664 



6.2. Finding background galaxies 



6.1. Observations and Data Reduction. 

NGC 3664 was observed with the WFPC2 in 1996 
March 3. The total exposure times are 600 seconds 
in the F555W passband, 800 seconds in the F814W 
passband, and 1600 seconds in the F439W passband. 
NGC 3664 is a Magellanic pecuhar (SBm(s)IV-V, de 
Vaucouleurs et al. 1991). From the recession veloci- 
ties (vgsr) cited in this reference for NGC 3664 and 
NGC 4536, and using the distance to NGC 4536 above 
(§13]), NGC 3664 is at a distance of 11.9 Mpc. The 
data will then have a spatial resolution of '--^ 6 pc 
per WF pixel. We assume the foreground Galactic 
absorption is As = 0.12 mag (de Vaucouleurs et al. 
1991). 

Figure |l^ shows the WFPC2 field, overlaid on an 
STScI Digitized Sky Survey image of the galaxy. Once 
again, we do not use the PC (WFl) frame. NGC 
3664 covers only 45% of the three WFC chips; this 
decreases the number of galaxies expected to be seen 
through the disk, but provides a fairly large surround- 
ing field for additional comparison. From its morphol- 
ogy, NGC 3664 is believed to be seen face-on. 

Partly because of the shallowness of the data, part- 
ly because the work on NGC 3664 started out in a 
relatively independent way, and also owing to a de- 
sire to explore diverse lines of attack, the analysis of 
the NGC 3664 data differs in the details from that of 
NGC 4536. The main differences are: (1) the subtrac- 
tion of the available (and properly scaled) i?F439W 
image from the / image to facilitate the search for 
relatively red objects (§|6.2|); (2) the use of the Rieke 
& Lebofsky (1985) Galactic reddening curve; (3) the 
semi-automatic identification of background galaxies, 
both real and simulated (|6.3|). To improve the statis- 
tics, we did more simulations at smaller extinction 
intervals (Table 3). We did not try to fix the warm 
pixels in the raw images, and cosmic rays were iter- 
atively removed based on the deviation of individual 
pixel values from the mean. Finally, we adjusted the 
parameters of SExtractor to minimize the scatter of 
the difference between the measured brightnesses of 
the HDF galaxies before and after adding them to 
the NGC 3664 fields. For the same reason, in what 
follows we use the "automatic" aperture magnitudes 
calculated by SExtractor, rather than the isophotal 



To remove most of the bright, blue star- forming 
regions in the disk of NGC 3664, we subtracted a 
properly-scaled F439W image from the F814W im- 
age. The scaling was done arbitrarily, to remove this 
blue component without creating substantial "holes" 
in the residual image. This method is reasonable, be- 
cause we cannot or can only barely see the background 
galaxies in the F439W passband. The disadvantage of 
this procedure is that it increases the noise. Howev- 
er, since background galaxies are identified primarily 
through their relatively red appearance, the disadvan- 
tage is far outweighed by removing the bright, rela- 
tively blue, foreground in the first step. The search for 
galaxies with SExtractor was performed in this arti- 
ficial, blue-subtracted image. The photometry of the 
identified galaxies, however, was done in the original 
F555W and F814W frames. 

6.3. The synthetic field simulations 

Table 3 shows the input parameters of the simu- 
lations, and the number and average output color of 
the HDF galaxies recovered behind and surrounding 
NGC 3664; a hue indicating the actual NGC 3664 
observations is included for comparison. As opposed 
to what we did in the case of NGC 4536, the selec- 
tion of galaxies out of the object catalogs produced 
by SExtractor was done in a semi-automatic fash- 
ion. For example, we only selected objects redder 
than {Vf555w — ^F814H') = 0.17 mag (in the cata- 
log, i.e., several tenths of mag redder in actuality in 
crowded regions, see §7.2). Other selection criteria 
were a (SExtractor) star-galaxy separation parameter 
"class_star" < 0.8 (the values go from 0.0 for an un- 
mistakable galaxy to 1.0 for a perfect star); a FWHM 
of a Gaussian fit to the object profile > 0"2; and 
an error in the "automatic" aperture magnitude < 
0.17 mag. In bright foreground star forming regions, 
the search and measurement of background galaxies is 
virtually impossible. Therefore, we excluded from the 
search regions of NGC 3664 with a surface brightness 
higher than 21 mag arcsec"^ (WF3) and 20.7 mag 



"Automatic" apertures are a constant, in this case 2.5, times 
the size of the Kron "radius", i.e., an elliptical aperture that 
contains the same fraction of the total light for all objects, 
regardless of their brightness profile. 
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arcsec"^ (WF4), respectively^. We still had to ex- 
clude by hand saturated stars that were included in 
the automatic selection of galaxies. 

6.4. Real vs. control galaxies 

Because of the difficulties in identifying galaxies 
reliably, we decided to exclude from the sample be- 
hind the disk of NGC 3664 all candidates fainter than 
'mpsiiw — 22.8. This reduced our sample of re- 
al galaxies found behind the projected disk of NGC 
3664 to 8 objects. To the same limiting magnitude of 
iTT-FsiiW = 22.8, in the simulation without extinction 
we recovered an average of 17 HDF galaxies within the 
disk. The discrepancy must be caused by extinction. 
For the areas outside the visible disk of the galaxy, 
we kept galaxies up to an integrated magnitude of 
mpsuw = 24. We found a total of 38 background 
galaxies in the surrounding field, vs. an average of 37 
control galaxies. 

6.5. Extinction in NGC 3664 

We find for the disk of NGC 3664 (Figure ||, top) 
an absorption of Aj ~ 1.01 ± 0.63 mag in the case 
of grey absorption, and Aj = 1.02 ± 0.63 mag in the 
case of a Galactic reddening law. The purely Poisson 
error is ±0.49 mag. For the surrounding area (Figure 
pjj , bottom) the results are, respectively, Aj < 0.31 
mag (< 0.17 mag without clustering error) for grey 
extinction and Aj < 0.30 mag (< 0.16 mag without 
clustering error) for Galactic reddening law. (The re- 
sults have been corrected for foreground extinction in 
the Milky Way.) Table 1 lists the relevant parameters 
for the calculation of the error in the extinction. 

6.6. Reddening in NGC 3664 

For the background galaxies, we measure a mean 
color of <{V - I)>inside= 1.46±0.19 behind the disk 
of NGC 3664, and of <(¥ - I)>outs^de^ 1-13 ± 0.10 
in the surrounding field. Figure ^ shows how we 
compare the color and extinction of the background 
galaxies to the simulations. The left panel exhibits 
the results for the field covered by the disk of NGC 
3664. 

As in the case of NGC 4536, the background galax- 
ies behind the disk of NGC 3664 lie in between the 



'These limits were set empirically, with the aid of the simula- 
tions, so as to retrieve the maximum number of HDF galaxies 
and the minimum amount of misidentifications 



lines for the grey and Galactic reddening simulations. 
To find the reddening, we compare the mean color of 
the background galaxies to the mean color of the grey 
simulation with the closest amount of extinction, i.e., 
the one with Ai = Ay = 0.95 mag. The output mean 
color of this simulation is <{V — I)>— 1.06 ± 0.11. 
Therefore, the background galaxies suffer a reddening 
of E{Vf555W - If814w) = 0.40 ± 0.22 mag. 

The right panel of Figure ^ shows the results 
for the galaxies recovered outside the visible disk of 
NGC 3664. The dashed line follows the color of the 
grey simulations after correcting for incompleteness: 
a shift was added to the average color of the recov- 
ered galaxies as a function of increasing extinction, so 
that it would follow the same behavior as the average 
color of the (un-added) HDF galaxies with a progres- 
sively brighter magnitude limit. Although we prefer 
to use the color without correction, because other- 
wise the small numbers make the error in the color 
very large, the line does indicate that we are seeing 
an incompleteness effect. Conversely, in the case of 
the simulations with Galactic reddening, the line after 
correction (not plotted) basically follows the output 
mean color of the recovered galaxies. This confirms 
that, given our color selection bias, we can still find 
relatively faint galaxies, even in the presence of ex- 
tinction, as long as they are reddened. The result 
for the background galaxies outside the visible disk of 
NGC 3664 is that their number and color are consis- 
tent with no attenuation. 

7. SYSTEMATIC ERRORS AND BIASES 

7.1. Color selection and the color magnitude 
relation 

Because spiral disks are relatively blue, red back- 
ground galaxies are more readily detected behind a 
foreground disk. As mentioned above, in crowded re- 
gions of NGC 4536 we selected in practice for red 
galaxies [i.e., with (Vpsssw — -^^^8141^) > 0.5], to be 
able to distinguish them from the stellar content and 
HII regions (which are predominantly bluer) in the 
foreground galaxy. This bias will be even more no- 
ticeable in grey simulations, where (given a certain 
attenuation at F555W) the F814W passband wiU be 
more attenuated than in the corresponding simulation 
with Galactic reddening. As a result, more faint blue 
galaxies will be discarded in the grey simulation than 
in the reddened one, the mean color measured for the 
outcome of the grey simulation will be redder than 
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otherwise, and the measured color shift between it 
and a reddened simulation will be smaller (bluer) than 
the input. The color-magnitude relation of galaxies 
(intrinsically brighter galaxies tend to be redder) will 
have an effect in the same direction: at higher extinc- 
tions, the statistically brighter galaxies will be detect- 
ed, and the average color of the sample will be redder. 

7.2. Systematic error in the color 

measurement due to the background 

As it searches for and measures objects, SExtrac- 
tor determines the background with a sigma-clipping 
algorithm, plus a modified mode in crowded areas. In 
very crowded fields, however, this results in an un- 
derestimation of the actual background when using 
large apertures (larger than 0"6 in radius for the two 
galaxies investigated here) . The background determi- 
nation algorithm has the purpose, precisely, of remov- 
ing bright spots (stars, clusters or Hll regions) from 
the measurement of the background, so that faint, ex- 
tended objects can be identified. On the other hand, 
the larger the photometric aperture, the higher the 
probability that these same bright spots in the fore- 
ground galaxy disk will lie within the aperture, and 
will be included in the measurement of the flux of the 
object of interest, in this case, the background galaxy. 

The difference between the true and the measured 
background is equal to the difference between the av- 
erage of the background and its mode. The result of 
the wrong background determination is equivalent to 
adding this difference to the galaxy luminosity mea- 
surement. We find that, on average, galaxies are 
measured 0.1 mag brighter in the F814W bandpass 
simulations when using isophotal magnitudes. In the 
F555W passband, the shifts depend on the crowd- 
ing. For NGC 4536, they go from ~ 0.1 in the 
WF4 frame to 0.6 mag for many galaxies fainter than 
mF555W = 22 in the WF2 frame, the most crowded of 
the three WFC frames. In short, fainter galaxies will 
be measured up to 0.5 mag too blue, if using isophotal 
apertures. 

As with the other errors discussed above, this is a 
systematic error that will affect both the simulations 
and the real galaxies, and therefore will not affect 
their mean relative reddening. However, for individ- 
ual galaxies, the fainter the galaxy the larger the er- 
ror will be (in magnitudes). Hence, this effect will 
increase the error in the estimate of the mean back- 
ground galaxy color. We have reduced this problem 
by determining the galaxy colors from a fixed aper- 



ture of 0''5 (5 pixel) diameter. The use of this aper- 
ture size decreased the shift between the input and 
the measured color in the simulations by a tenth of a 
magnitude, and the dispersion in the shift by up to 
a tenth of a magnitude as well. The error introduced 
by not measuring possible color gradients over larger 
background galaxy radii should be much smaller than 
this error, especially when the measurement of a color 
gradient depends so critically on an accurate estima- 
tion of the background. Much more importantly, by 
choosing an aperture to perform the color measure- 
ment, we are choosing the scale over which w e ar e 
averaging the effects of the obscuring material (§8.3). 
At the distance of NGC 4536, 0^5 is ~ 40 pc; for 
NGC 3664, the hnear size of the aperture is ^ 30 pc. 

7.3. Systematic errors and the derived 
extinction law 

The systematic errors discussed above in §^!J and 



7.2 will affect the simulations and the real galaxies 



alike, and hence will not have an effect on the redden- 
ing determination. The only error that could have an 
impact on the determination of the reddening is the 
interpolation error in the HDF Vfsssw images (§2.1). 
However, the mean expected shift of ~ -0.03 mag be- 
tween the interpolated and the real tofsssvk of the 
HDF galaxies is smaller than (1) the dispersion of the 
shift; (2) the error in the mean colors of the real galax- 
ies and of the sets of output simulated galaxies (§7.5); 
and (3) the dispersion in the measurements of individ- 
ual galaxy colors(§5.2). We conclude that the error 
introduced by the interpolation will have a negligible 
effect over both the background galaxy detection and 
the reddening determination. 

7.4. The clustering error 

To estimate the contribution from galaxy cluster- 
ing to the error in the number-counts, we take the 
two-point correlation function w{9) from Roche et al. 
(1993) or from Brainerd et al. (1995), depending on 
the completeness magnitude limit of the data. We 
assume (i? — /) ~ 0.5 (Small et al. 1995), since these 
references list w{6) for limits in the i?-band. Table 1 
shows, for each disk (NGC 4536 and NGC 3664) and 
field, the angular size 9 of the field, the assumed (ap- 
parent) R magnitude completeness limit, w{9)^ the 
number of background galaxies, the derived Ai (red- 
dened and grey) extinction, the total error in the ex- 
tinction, the error contributed by the Poisson error in 
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the number-counts, and the error owing to clustering 
errors. 

To estimate how the simulations contribute to the 
clustering error, rather than presenting data for each 
of the two simulations used in the calculation of the 
extinction (the simulations with numbers of recovered 
galaxies immediately above and below the number of 
real galaxies), we list data for only one fiducial simu- 
lation with the same number (per WFC field) of re- 
covered HDF galaxies as the real background galax- 
ies. Even though the 3 HDF WFC fields have been 
added one at a time to each of the WFC fields of 
the investigated galaxies, in reality they are of course 
contiguous, and have to be treated as one single field 
for the calculation of the clustering error. The errors 
contributed by the simulations to the final estimation 
of the extinction will be \/2 the numbers given for the 
fiducial simulation. Except for NGC 4536, WF4, all 
the other values of w{d) (and therefore of the cluster- 
ing errors) are slightly overestimated, both because 
we take the correlation function for a square with 
an area equal to the area under consideration (not 
a square), and because Roche et al. (1993) include 
galaxies with 18 < rrifl < 20 in their calculation of 
the correlation function; these galaxies add very sig- 
nificantly to wlO), but there are none in our sample of 
background objects (we detect no more than a couple 
in the simulations without extinction) . Conversely, in 
the case of the NGC 4536 WF2 and WF3 fields (an- 
alyzed in a combined fashion), there is a correlated 
error in the simulations that we have not included, 
and that is due to the fact that we count the same 
HDF fields twice. 

7.5. Photometric error 

If the color distribution of galaxies was Gaussian 
and we were sampling it in a systematic way, we could 
confidently equate the error in the mean color to its 
measured dispersion divided by \/iV, where N is the 
number of galaxies used for the color measurement. 
However, neither of these assumptions is true. The 
color distribution of galaxies depends on their lim- 
iting magnitude and has an intrinsic dispersion that 
is not Gaussian. Moreover, in our case we have to 
infer their mean color from the few galaxies that we 
can find in the presence of crowding. Crowding will 
have three adverse effects on the color measurement. 
The color-magnitude distribution of galaxies resem- 
bles an inverted triangle (cf. for example Williams et 
al. [1996]). (1) Crowding will push the completeness 



limit up to brighter magnitudes, where the intrinsic 
color dispersion of galaxies is largest. (2) It will re- 
duce the number of galaxies available for the analysis 
in a dramatic way (Figure ^). (3) Crowding will de- 
plete galaxies in an inhomogeneous fashion: in combi- 
nation with clustering, it will determine the random 
sample of background galaxies from which we will try 
to infer the mean color of all background galaxies. 

To check our estimate of the error in the color 
measurement, we have looked at the dispersion in 
the mean color of 100 random samples of N galax- 
ies (where N is equal to the number of background 
galaxies) taken from the HDF galaxies, to our mag- 
nitude limit. We have repeated the procedure with 
the simulations, and we have found that the predict- 
ed error is of the same order as the error found by 
assuming a Gaussian distribution. Most of the time, 
the error from bootstrap statistics (Babu & Feigel- 
son 1996) is slightly smaller than the Gaussian error, 
probably because the latter is calculated from the real 
sample and therefore includes the effect of clustering 
and other systematics. In the case of NGC 4536, we 
adopt the larger of the two error estimates for each 
simulation, as well as for the real galaxies. In the case 
of NGC 3664, we use the Gaussian uncertainty. 

8. SUMMARY AND DISCUSSION 
8.1. NGC 4536 

Summarizing the results (errors in parenthesis are 
Poisson only, no clustering) : the northwestern arm of 
NGC 4536 shows an attenuation of Ai = 0.74 ± 0.49 
(± 0.26) mag (grey extinction) or Aj = 1.07 ±0.49 
(± 0.26) mag (Galactic reddening law), and a red- 
dening of E{Vf555W - Ifsi4w) < 0.32 mag (grey ex- 
tinction) or E{Vf555w — If814w) < 0.29 mag (Galac- 
tic reddening law). The +1<t error of 0.18 mag is 
already included in the color uncertainty. For the 
northern interarm and outside regions, the results are, 
respectively, Aj < 0.46 (< 0.33) mag (grey extinc- 
tion) or A] < 0.52 (< 0.39) (Galactic reddening), and 
E{Vf555w — ^F8i4H') = 0.20 ± 0.13 mag. The extinc- 
tion values for the interarm region already include a 
la error of -1-0.37 (0.21) mag. The Poisson errors in 
parenthesis would be adequate for the extinction if 
the galaxy density was the same behind NGC 4536 as 
it is in the HDF. 



14 



8.2. NGC 3664 

To summarize the results: the disk of NGC 3664 
shows an attenuation of Ai = 1.02 ± 0.63 (± 0.49) 
mag (Galactic reddening law) or Aj = 1.01 ± 0.63 
(± 0.49) mag (grey extinction), and a reddening of 
E{Vf555W - If814w) = 0.40 ± 0.22 mag. The Pois- 
son errors in parenthesis would be adequate for the 
extinction if the galaxy density was the same behind 
NGC 3664 as it is in the HDF. As it turns out, in the 
surrounding field of NGC 3664, this is precisely the 
case. 

8.3. Conclusions and discussion 

We have used background galaxy counts and col- 
ors, for the first time, to measure the extinction and 
reddening through two very different galaxy disks. 
The images also have very different available expo- 
sure times, which allowed us to test the limits of the 
method. 

Our principal conclusions are listed below. 

(1) In the case of spiral disks, the background 

galaxy approach is mostly limited by crowding and 
confusion in the foreground disk, even more than by 
the background galaxy clustering. We have devel- 
oped and calibrated the "synthetic field" method to 
decouple the effects of crowding from those of extinc- 
tion. The method comprises the production of arti- 
ficial frames by adding suitable and properly scaled 
reference fields, in this case the HDF, directly into 
the spiral galaxy images. The reference frames are 
attenuated to mimic different amounts of extinction 
and reddening; the idea is to attenuate the control 
field until one recovers the same number of simulat- 
ed galaxies as there arc real backgroimd galaxies. To 
ensure that systematics and confusion affect synthet- 
ic and real galaxies in the same way, both real and 
control galaxies are identified in the same fashion 
(i.e., the known positions of HDF galaxies are nev- 
er used). The attenuation derived from the number- 
counts should correspond to the average extinction 
through the disk of the galaxy. Finally, a comparison 
between the average colors of, respectively, the con- 
trol HDF galaxies and the real background galaxies 
will allow an estimate of the degree of reddening. 

(2) We tried slightly different lines of attack with 
both galaxies, but found them to be fairly equivalent. 
When a B image is available, like for NGC 3664, its 
subtraction from the I image does seem to facilitate 
the search for galaxies, at least in frames with short 



exposiire times. It would also be desireable to auto- 
mate as much as possible the final identification of 
real and synthetic galaxies. Unfortunately, a fully 
automatic identification might be hard to achieve in 
frames as crowded as the ones we are studying. 

(3) There is a significant amount of dust in NGC 
3664, even though it is a Magellanic irregular. 

(4) There is more extinction in the northwestern 
arm of NGC 4536 than in the northern interarm re- 
gion of the galaxy. The reddenings, however, are com- 
parable. We have mentioned before that we derive 
both parameters independently: the extinction from 
the galaxy number-counts, and the reddening from 
the background galaxy photometry. The method we 
use to measure the extinction will take into account 
galaxies that are so attenuated that we cannot see 
them, while we measure the reddening from galax- 
ies that are attenuated only to a point where we can 
still observe them. Also, although we average the col- 
or excess over a number of background galaxies, the 
color of each one of these galaxies is measured in an 
aperture with a linear size of ^ 40 pc at the distance 
of NGC 4536. It is possible that at this scale the 
dust distribution is genuinely smooth (Keel & White 
1997). 

In view of all the above, in principle we prefer to 
present the results obtained from the number-counts 
and from the colors independently, and not in the 
form of a reddening law. The number-counts measure 
the total average extinction through a whole region; 
the colors measure the extinction owing to a diffuse 
component. 

On the other hand, if we assume that nowhere in 
the region of interest the extinction is so high that 
we cannot sample it with background galaxies, it is 
possible to determine an extinction law from the two 
values (Appendix). This seems to be a safe assump- 
tion for the northern interarm region of NGC 4536, 
and there both the extinction and the reddening are 
consistent with a Galactic reddening law. The sim- 
ilar reddening but significantly higher extinction in 
the northwestern arm can be explained by the com- 
bination of a diffuse obscuring component with re- 
gions of higher opacity. If, in spite of the caveats, we 
use the extinction and the reddening measured in the 
northwestern arm to derive a reddening law, then the 
extinction there is greyer than Galactic and can be 
caused by unresolved clumps (Witt & Gordon 1996), 
in this case at the scale of ~ 40 pc. 
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(5) For the interarm region of NGC 4536, the ex- 
tinction can be deprojected to a face-on Aj < 0.22 
mag. For the northerm arm of NGC 4536, however, 
given the patchy dust distribution, it is not possi- 
ble to derive a number for the face-on extinction in 
a model-independent way. We can, however, assume 
that the reddening of the background galaxies behind 
the arm is caused by a smooth component and follows 
a Galactic extinction law. In that case, the redden- 
ing implies a face-on extinction through the diffuse 
component in the arm of < 0.16 mag. 

(6) The usefulness of the background galaxy counts 
and colors to measure foreground extinction has been 
recently confirmed by their application to the field 
of the gamma-ray burster GRB970228 (Fruchter et 
al. 1998). In that analysis, we used four different 
reference frames to measure the extinction in the 
F606W passhand: the HDF; the weak radio galaxy 
53W002 field; and two fields located, respective- 
ly, at J(2000) 15:58:49.8, 42:05:23, and at J(2000) 
14:17:43.63, 52:28:41.2 (Westphal, Kristian, & Groth 
1994). Since the GRB970228 field is not crowded, 
it was not necessary to produce "synthetic fields". 
Under those circumstances, and assuming a Galactic 
reddening law, both counts and colors yield the same 
result (respectively, 0.42±0.23 mag, and 0.54±0.04 
mag) . The average of these two measurements of the 
Galactic extinction towards the gamma-ray burster 
(0.48±0.12 mag) is consistent with the value of 0.72 
±0.24 mag obtained by Burstein & Heiles (1982), and 
with the extinction of 0.69±0.11 measured from maps 
of dust IR emission (Schlegel, Finkbeiner, & Davis 
1997) . The marginal systematic difference in the mea- 
surements can be explained by our better spatial res- 
olution (roughly 1', vs. 6!l of Schlegel et al. and 0?6 
of Burstein & Heiles), since the burster is located at 
a position {I = 188.913°, b = -17.941°) where the 
Milky Way exhibits a steep gradient in extinction (cf. 
Burstein & Heiles 1982). 

(7) From the Cepheids detected in the same NGC 
4536 data, Saha et al. (1994) derive a mean extinction 
of 0.07±0.05 mag in Kron-Cousins Likewise, they 
estimate the / differential extinction between the arm 
and the interarm regions to be 0.05±0.10 mag. These 
numbers imply for the arm of NGC 4536 an extinc- 
tion ^7^0.2 mag, much lower than what we find from 



From Figure 12 of Holtzman ct al. (1995), the extinction in 
Kron-Cousins / is directly comparable with the extinction in 



the galaxy number-counts. More importantly, the de- 
tected Cepheids span a range of at least 2 magnitudes 
in brightness; Saha and co-workers had the capability 
to measure an extinction much higher than 0.2 mag. 
The fact that they did not implies that the Cepheids 
detected in the arm lie preferentially either in areas 
of low extinction or on top of the obscuring material. 
Together with the extinction derived from the galaxy 
number-counts, the Cepheid reddening is consistent 
with the existence of regions of significant opacity in 
the arm of NGC 4536. On the other hand, the extinc- 
tion measured from the Cepheids agrees with the one 
we find for the interarm region and, from the galaxy 
colors, for the diffuse component in the arm (<0.52 
and <0.42 mag, respectively, for a Galactic reddening 
model). Lastly, since the Cepheids are embedded in 
the disk, statistically they should suffer only about 
half of the total extinction, which background illu- 
minating objects can probe; our measurements and 
those by Saha and collaborators are consistent with 
this picture. 

(8) That in the case of spiral galaxy disks we are 
totally dominated by crowding is demonstrated by the 
fact that there is less than 30 percent improvement in 
the extinction errors for NGC 4536, relative to NGC 
3664, even though we have 5 times better signal-to- 
noise ratio in Ifbiaw and more than 10 times better 
signal-to-noise ratio in VpbbbW- Likewise, the color 
errors are only 20 percent better for NGC 4536. Per- 
haps even more telling, the accuracy of both the ex- 
tinction and the color measurements is the same in 
the less obscured regions of the galaxies, despite the 
enormous difference in exposure times. 

We can both reduce the contribution to the rela- 
tive errors from the control fields, and perhaps even 
measure the clustering error directly, by using an in- 
creasing number of them. And since the dominating 
source of error is crowding, it certainly will not mat- 
ter if we use reference fields with shorter exposure 
times than the HDF to produce the "synthetic fields" . 
However, even if we manage to determine very accu- 
rately the mean and the dispersion of the number of 
galaxies in a typical WFPC2 reference field, we will 
still be left with the larger source of error: the stud- 
ied galaxies themselves, with their "single realization" 
background fields. If NGC 4536 and NGC 3664 are 
representative, the extinction errors will go from 0.3 
mag in galaxies with exposures of several x 10* sec- 
onds, to 0.5 mag in galaxies with shallow exposures of 
a few X 10'^ seconds; the error in the color will amount 
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to ^ 0.1 mag in all cases. In the future we will try 
to reduce this source of error by combining together 
the results from background galaxies behind spirals of 
similar Hubble type, inclination, and distance to the 
Milky Way. 
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APPENDIX 
METHOD-DEPENDENT GREYNESS? 
The reddening and the attenuation are derived independently with our approach. In what fohows, we discuss 
the hmitations of the method for the derivation of a meaningful extinction law. 

Let us consider the case where each of the background galaxies used to derive the "average" extinction through 
the disk suffers from a different amount of Galactic reddening. Our method will derive a reddening law that can 
be described by the following equation: 

17= , (1) 

E{y~i) ^' 



where Ai denotes the extinction derived from the number-counts and E{V — /) is the average reddening of the 
recovered galaxies. E{V — I) can also be written as: 

i i 

where the summation encompasses all the background galaxies. 
The extinction, on the other hand, can be expressed as i^^- 

N 

Aj = ~-2.5-C-logi — ), (3) 

where C describes selection biases and A^'o is the normalization. 

If P{mi, All) is the probability (0 or 1) of finding a galaxy with magnitude rrii and suffering an attenuation An, 
we can rewrite the above equation as: 

Ao Y 

To obtain the attenuation in a fractional area fj, with attenuation Ajj, of the original image, the equation 
becomes: 

Ai, = -2.5 • C ■ /og(-i- 5^P(m„ A,,)). (5) 

Splitting the summation in equation (4) into two summations, one running over the areas with a given attenuation 
Aj, the other running over all galaxies in each of these areas, equation (4) becomes: 

17=-2.5.C-Mi- ^ (6) 

area.j i 

Using the value of P{mi, ^^^j) from equation (5), and assuming that J2area j fj ~ finahy obtain: 

A;^-2.5-C-log{Y, fA^"^^'')- (7) 

area,j 

The "reddening law" would read: 



Ih = — ^ -Ri (8) 
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Equation (7) shows that the extinction obtained from the number-counts is a hnear average over the whole 
area considere(^, and as such takes into account real background galaxies that we cannot see, while the color 
excess involves the average of the absorption in magnitudes, and only over background galaxies that we actually 
see. Strictly, the measured Rj and the intrinsic Rj will be equal only if Aj is the same for all galaxies. But, 
mainly, it is possible to measure a much higher extinction from the number-counts than the extinction implied by 
the reddening. 



That is, the whole area that is sampled by the reference galaxies, so crowding will not affect the extinction measurement 
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Table 1 
Attenuation Results. 



Object/Field 


e 


R limit 


w{e) 


N 


Al,Gal 


-^I,grey 




(^Aj,Poisson 




NGC 4536, WF2 & 3 


113" 


22.2 


0.0231^ 


29 


1.07 


0.74 


0.49^ 




0.20 




0.25 


Fiducial Simulation 


139" 


22.2 


0.020'" 


87 










0.12 


V2x 


0.23 


NGC 4536, WF4 


80" 


24.2 


0.012^= 


35 


<0.52^ 


<0.46^ 


0.37*' 




0.18 




0.18 


Fiducial Simulation 


139" 


24.2 


0.008<= 


105 








V2x 


0.11 


V2x 


0.15 


NGC 3664, inside^ 


95" 


22.3 


0.026'" 


8 


1.02 


1.01 


0.63^ 




0.38 




0.26 


Fiducial Simulation 


165" 


22.3 


0.016'" 


24 








V2x 


0.22 


V2x 


0.21 


NGC 3664, outside^ 


101" 


23.8 


0.014'^ 


38 


<0.30'= 


<0.31*' 


0.37^ 




0.18 




0.19 


Fiducial Simulation 


175" 


23.8 


0.009'^ 


114 








x/2x 


0.10 


x/2x 


0.15 



Cols.- (1) Object and field; (2) angular size of the field; (3) (unextincted) R completeness limit (apparent 
mag), assumed to be 0.5 mag fainter than the mF8i4,w limit; (4) two-point correlation function; (5) number of 
background galaxies; (6) Aj (mag) inferred from N, assuming a Galactic reddening law; (7) Aj assuming grey 
extinction; (8) total error in Aj (mag); (9) error from the shot-noise in the number counts; (10) error contributed 
by galaxy clustering. 

^Aj values include a correction for Galactic extinction of Aj^MUkyWay — 

0.04, extrapolated from = 0.12 (de 

Vaucouleurs et al. 1991; Rieke & Lebofsky 1985). 
'"Roche et al. 1993. 
"^Brainerd et al. 1995. 

"^Interpolated between values in the two previous references. 

''Value already includes +la. 

^Error includes contribution from simulations. 
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Table 2 
Simulations for NGC 4536. 



Num/Type 


Ai 


Av 






<V -I >2&3 


<V -I>i 


1, no attenuation 


0.00 


0.00 


168 ± 38 


113 ± 18 


1.11 ± 0.07 


1.20 ± 0.06 


2, grey 


0.58 


0.58 


99 ± 23 


60 ± 11 


1.20 ± 0.06 


1.17 ± 0.07 


3, mixed 


0.58 


0.79 


111 ± 26 


65 ± 12 


1.36 ± 0.05 


1.36 ± 0.06 


4, Galactic reddening 


0.58 


1.02 


117 ± 27 


81 ± 14 


1.53 ± 0.07 


1.60 ± 0.07 


5, grey 


1.02 


1.02 


65 ± 16 


43 ± 9 


1.33 ± 0.07 


1.25 ± 0.10 


6, grey 


1.51 


1.51 


36 ± 10 


27 ± 6 


1.29 ± 0.09 


1.37 ± 0.12 


7, mixed 


1.51 


2.06 


49 ± 13 


36 ± 8 


1.78 ± 0.08 


1.83 ± 0.10 


8, Galactic reddening 


1.51 


2.60 


60 ± 15 


47 ± 9 


2.10 ± 0.10 


2.04 ± 0.12 


9, grey 


2.60 


2.60 


14 ± 5 


12 ± 4 


1.35 ± 0.18 


1.45 ± 0.23 


NGC 4536 






87 ± 26^ 


105 ± 25^ 


1.36 ± 0.14 


1.40 ± 0.12 



Cols.- (1) Simulation number and type; (2) input attenuation (mag) in the F814VF passband; (3) input at- 
tenuation (mag) in the FbbbW filter; (4) recovered galaxies in the WF2 and WF3 fields (errors include Poisson 
and clustering uncertainty, see §7.4 and Table 1); (5) recovered galaxies in the WF4 field; (6) mean color of the 
recovered galaxies, WF2 and WF3; (7) mean color of the galaxies recovered in the WF4. 

^The number-counts have been multiplied by 3 for easy comparison with the simulations. 
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Table 3 
Simulations for NGC 3664. 





Nuni/T\"])(> 


Ai 




71 






A 






< V - I >,„ 


< V - I >„„t 


1, 


no attenuation 


0.00 


0.00 


50 


± 


12 


112 


± 


19 


0.97 ± 0.07 


1.08 ± 0.05 


2, 


grey 


0.15 


0.15 


47 




11 


92 


± 


16 


0.98 ± 0.08 


1.12 ± 0.06 


3, 


Galactic reddening 


0.15 


0.31 


48 


± 


11 


103 


d= 


18 


1.13 ± 0.07 


1.27 ± 0.07 


4, 


grey 


0.40 


0.40 


39 


± 


10 


68 


± 


13 


1.00 ± 0.09 


1.21 ± 0.07 


5, 


Galactic reddening 


0.40 


0.83 


42 


± 


10 


86 




15 


1.57 ± 0.09 


1.47 ± 0.07 


6, 


grey 


0.75 


0.75 


33 


± 


8 


53 


± 


10 


1.01 ± 0.09 


1.23 ± 0.07 


7, 


Galactic reddening 


0.75 


1.56 


36 


± 


9 


62 




12 


1.69 ± 0.10 


1.79 ± 0.09 


8, 


grey 


0.95 


0.95 


28 


± 


7 


42 


± 


9 


1.06 ± 0.11 


1.22 ± 0.09 


9, 


Galactic reddening 


0.95 


1.97 


28 


± 


7 


46 


± 


9 


1.82 ± 0.12 


1.86 ± 0.11 


10, 


grey 


1.15 


1.15 


20 


± 


6 


43 


± 


9 


1.13 ± 0.13 


1.20 ± 0.09 


11, 


Galactic reddening 


1.15 


2.39 


21 


± 


6 


37 


± 


8 


1.99 ± 0.14 


2.12 ± 0.08 


12, 


grey 


1.30 


1.30 


12 


± 


4 


39 


± 


8 


1.01 ± 0.17 


1.21 ± 0.10 


13, 


Galactic reddening 


1.30 


2.70 


17 


± 


5 


33 


d= 


7 


2.12 ± 0.17 


2.36 ± 0.11 




XGC 366 1 






21 


± 


10" 


111 


± 


27a 


I.IG ± 0.19 


1.13 ±0.10 



Cols.- (1) Simulation number and type; (2) input attenuation (mag) in the F814W passband; (3) input atten- 
uation (mag) in the F555W filter; (4) recovered galaxies within the disk of NGC 3664 (errors include Poisson and 
clustering uncertainty see §7.4 and Table 1); (5) recovered galaxies surrounding the disk of NGC 3664; (6) mean 
color of the recovered galaxies within the disk; (7) mean color of the galaxies recovered outside. 

^The number-counts have been multiplied by 3 for easy comparison with the simulations. 
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Fig. 2. — Cumulative histograms of galaxy counts in Ipsuw images. Left: WF2 and WF3. A: Field galaxies in 
the "degraded" HDF; b: HDF galaxies recovered from the synthetic fields without extinction; c; HDF galaxies 
recovered from the simulation with 0.6 mag of extinction (Galactic reddening curve) at IpsiiW] d: Real background 
galaxies in the NGC 4536 frames. The numbers of galaxies have been normalized to the area of one WFC field. 
Right: WF4. Histograms labeled as in left panel (no simulation with extinction is included). 
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Fig. 3. Number of HDF galaxies per (isophotal) mpsnw bin. The HDF's signal-to-noise ratio has been degraded 
to the level adequate to the exposure time and sky background of the NGC 4536 image. Top left: WF2; top right: 
WFS; bottom left: WF4; bottom right: average of the 3 WF fields. In general, comparisons of real and simulated 
background galaxies will be done against this average. 
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Fig. 4. — Same as figure for the Vfsssw filter. The galaxies were searched for in the Ipsuw images, so they are 
the same objects shown in figure 13, and their light has been integrated up to the same limit. 
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Fig. 5. — STScI Digitized Sky Survey image of NGC 4536, displayed with a linear scale. The blue (J) Science and 
Engineering Research Council Survey plate was obtained on 1979 May 19, with the UK Schmidt Telescope. The 
position of the WFPC2 field is superimposed. North is up and East is to the left. 
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Fig. 6. — Cumulative histograms of real background galaxies found in the NGC 4536 / images. Left: WF2 and 
WF3. The numbers have been normalized to the area of one WFC field. Right: WF4. 
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Fig. 7. — Cumulative histogram of the galaxies in the "degraded" HDF Ifsuw images. Top left: WF2; top right: 
WF3; bottom left: WF4; bottom right: average of the 3 WF fields. The differences in the number-counts owing to 
galaxy clustering are most noticeable at (isophotal) magnitudes brighter than mpsMW ~ 24. 
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Fig. 8. — Comparison between the average number of HDF galaxies recovered from the NGC 4536 frames in the 
synthetic field without extinction (solid line), and the average number per (isophotal) mpsuw bin of all HDF 
galaxies (dashed line). The HDF's signal-to-noise ratio has been degraded to the level adequate to the exposure 
time and sky background of the NGC 4536 image. The figure shows the dramatic effect of crowding alone on the 
number counts. Left: WF2 and WF3. Given the isophotal brightness limit, while the HDF galaxies are complete 
to at least mpsuw = 25, when added to the NGC 4536 they are complete only to mpsuw = 23. The numbers 
have been normalized to the area of only one WF field. Right: WF4. Simulation is complete to rupsuw = 24. 
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Fig. 9. — Same as figure for the Vpsssiv passband. The objects in the synthetic fields were searched for in the 
If814 frames, so they are the same as those shown in figure ^ Many more faint, blue galaxies were detected in the 
WF4 NGC 4536 frame (nght), relative to the WF2 and WF3 fields combined {left). 
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Fig. 10. Comparison between the number of baekground galaxies {solid line) and the average number per 
(isophotal) mF8i4w bin of HDF galaxies recovered from the NGC 4536 frames in the synthetic field without 
attenuation (dashed line). Left: WF2 and WF3. The effects of extinction are readily noticeable, both from the 
reduction in the number-counts and from the shift to fainter magnitudes of the real galaxies. The numbers have 
been normalized to the area of only one WF field. Right: WF4. Any effects are much less apparent. 
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Fig. 11. — Comparison between the number of background galaxies {solid line) and the average number per mF555W 
bin of HDF galaxies recovered from the NGC 4536 frames in the synthetic field without attenuation {dashed line). 
The galaxies were searched for in the If814W images, and hence they are the same as those in figure |l^. Left: 
WF2 and WF3. Both the numbers and the brightnesses of the background galaxies are reduced, compared to 
the simulation. The numbers have been normalized to the area of only one WF field. Right: WF4- Effects of 
extinction, especially in the numbers, are much less apparent. 
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Fig. 12. — Comparison between tlie (Vpsssiy — If814w) color of the bacliground galaxies {solid line) and the color 
of the HDF galaxies recovered from the NGC 4536 frames in the synthetic field without attenuation {dashed line). 
Left: WF2 and WF3. Their numbers look depleted, but the effects of extinction on the color of the real galaxies 
are not very noticeable. The numbers have been normalized to the area of only one WF field. Right: WF4.. Both 
the numbers and the colors of the real galaxies are consistent with those of the simulation. 
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Fig. 13. — Top: NGC 4.536, WF2 and WF3 frames. Filled triangles and solid line: average number of recovered 
HDF galaxies, vs. simulated Galactic reddening extinction of Ai. Open triangles and dotted line: average number 
of recovered HDF galaxies vs. extinction, when using a grey extinction model for the simulations. The difference 
between the extinction inferred from the two models increases with extinction. The error bars include the statistical 
Poisson uncertainties, as well as the field-to-field variations expected from galaxy clustering. Horizontal solid line: 
number of real background galaxies with its uncertainty [horizontal heavy dashed lines). Bottom: NGC 4536, WF4 
frame. Symbols as in top panel. 
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Fig. 14. — Left: results for NGC 4536, WF2 and WF3 fields. Open triangles and dotted line: dependence of 
the mean measured color of the simulations on the applied grey absorption. Filled triangles and dotted line: mean 
measured color of the simulations with extinction following a Galactic reddening law. Solid line: theoretical (input) 
Galactic reddening law. Open square: location of the background galaxies behind NGC 4536, derived from a grey 
absorbed model. Filled square: location of real background galaxies behind NGC 4536, using a Galactic reddening 
law. Right: Results for NGC 4536, WF4 field. Symbols as in left panel. 
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Fig. 15. — Mixed simulations for NGC 4-536. Left: WF2 and WF3. Open circles and thick-dotted line: mixed 
simulations. Thick-solid line: input of the mixed simulations. Thin-solid lines: results of Galactic reddening and 
grey simulations (Figure [l^). Right: WF4. Symbols as in left panel. 
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Fig. 16. — STScI Digitized Sky Survey image of NGC 3664 displayed with a logarithmic scale. The red (E) National 
Geographic Society-Palomar Observatory Sky Survey plate was obtained on 1955 April 16, with the 48-inch Oschin 
Schmidt Telescope. The position of the WFPC2 field is superimposed. North is up and East is to the left. 
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Fig. 17. — Top: results for the disk of NGC 3664- Open triangles and dotted line: average number of recovered 
HDF galaxies vs. simulated grey extinction of Aj. Full triangles and solid line: average number of recovered HDF 
galaxies vs. extinction, when using a Galactic reddening model for the simulations. The error bars include Poisson 
uncertainty in the number-counts and clustering error. Horizontal solid line: number of real background galaxies 
with its Poisson and clustering uncertainty {horizontal short-dashed lines). Bottom: results for the surrounding 
field. Symbols as in top panel. 
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Fig. 18. — Left: results for area inside the disk of NGC 3664- Open triangles and dotted line: dependence of the 
mean measured color on the applied extinction for the simulations with grey absorption. Filled triangles and dotted 
line: mean measured color of the simulations with Galactic reddening. Solid line: theoretical (input) Galactic 
reddening law. Open square: location of the background galaxies inside the disk area of NGC 3664. Right: results 
for area outside of the disk of NGC 3664- Long-dashed line: measured colors of the grey simulations, corrected for 
incompleteness. Symbols as in left panel. 
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